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ABSTRACT 

The Sersic model has become the standard to parametrize the surface brightness distribution of early-type galaxies and bulges of spiral galaxies. 
A major problem is that the deprojection of the Sersic surface brightness profile to a luminosity density cannot be executed analytically for 
general values of the Sersic index. Mazure & Capelato (2002) used the Mathematica® computer package to derive an expression of the Sersic 
luminosity density in terms of the Meijer G function for integer values of the Sersic index. We generalize this work using analytical means and 
use Mellin integral transforms to derive an exact, analytical expression for the luminosity density in terms of the Fox H function for all values of 
the Sersic index. We derive simplified expressions for the luminosity density, cumulative luminosity and gravitational potential in terms of the 
Meijer G function for all rational values of the Sersic index and we investigate their asymptotic behaviour at small and large radii. As 
implementations of the Meijer G function are nowadays available both in symbolic computer algebra packages and as high-performance 
computing code, our results open up the possibility to calculate the density of the Sersic models to arbitrary precision. 
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1. Introduction 

The Sersic model (Sersic 1968) is a three-parameter model for 
the surface brightness profile of galaxies that has been intro- 
duced as a generalization of the de Vaucouleurs R 1 ^ 4 model 
(de Vaucouleurs 1948). It is defined as 



I(R) = 7 exp 



l/m 



(1) 



where I(R) is the surface brightness at radius (on the plane of 
the sky) R, Iq is the central surface brightness, R e is the effec- 
tive radius, and m is the so-called Sersic index that describes the 
index of the logarithmic slope power law. The parameter b is 
a dimensionless parameter that depends on the Sersic index m 
and whose value can be derived from the requirement that the 
isophote corresponding to R e encloses half of the total flux. 

Over the past two decades, the Sersic model has be- 
come the standard to describe the surface brightness pro- 
files of early-type galaxies and bulges of spiral galaxies (e.g. 
Daviesetal. 1988; Caon et al. 1993; D'Onofrio et al. 1994; 
Cellone et al. 1994; Andredakis et al. 1995; Prugniel & Simien 
1997; Mollenhoff & Heidt 2001; Graham & Guzman 2003; 
Allen et al. 2006; Gadotti 2009). In the past few years, Sersic - 
like models have also gained popularity as a model to describe 
the spherically averaged profiles for dark matter haloes. While 
models with a power law behaviour at small and large radii 
were preferred in earlier simulations (e.g. Navarro et al. 1997; 
Moore et al. 1999), other models seem to fit the mass den- 
sity distribution of more recent (and higher resolution) sim- 
ulations better (Navarro et al. 2004, 2010; Merritt et al. 2005, 
2006; Graham et al. 2006; Aceves et al. 2006; Gao et al. 2008; 
Duffy et al. 2008; Stadel et al. 2009). Among these models, also 
the Sersic model (i.e. a model where the projected surface den- 
sity is described by a Sersic law) has been proposed as a univer- 
sal description for simulated dark matter haloes. 

Mainly as a result of its popularity to describe the sur- 
face brightness profiles of early-type galaxies, the properties 



of the Sersic model have been examined in great detail. Ciotti 
(1991) and Ciotti & Lanzoni (1997) give a detailed descrip- 
tion of the properties of the Sersic model, including spatial 
and dynamical properties. Ciotti & Bertin (1999) give a full 
asymptotic expansion of the dimensionless scale factor b of the 
Sersic model. Graham & Driver (2005) present a compendium 
of mathematical formulae on photometric parameters such as 
Kron magnitues and Petrosian indices, and Cardone (2004) and 
EKasdottir & Moller (2007) investigate the lensing properties. A 
major problem with the Sersic models is that the deprojection of 
the surface brightness profile to a luminosity density is in gen- 
eral non-analytical. In practice, one often uses approximations 
for the Sersic models when the luminosity density (or the mass 
density when the Sersic model is used to describe the distribu- 
tion of dark matter) is necessary (e.g. Prugniel & Simien 1997; 
Lima Neto et al. 1999; Trujillo et al. 2002). An unexpected an- 
alytical progress was the work by Mazure & Capelato (2002), 
who demonstrated that it is possible to elegantly write the spa- 
tial luminosity density of the Sersic model in terms of the Meijer 
G function. Unfortunately, their result only holds for integer val- 
ues of the Sersic index, which is a significant limitation for prac- 
tical applications. Moreover, since their result fell as a deus ex 
machina out of the Mathematica® computer algebra package, it 
is hard to see whether it can be generalized to all Sersic indices. 

In this paper, we tackle the deprojection of the Sersic surface 
brightness profile using analytical means. We apply an integra- 
tion method based on Mellin integral transforms and derive an 
exact, analytical expression for the luminosity density in terms 
of the Fox H function for all values of the Sersic index m. For ra- 
tional values of m, the luminosity density can be written in terms 
of the Meijer G function. As the Meijer G function is nowadays 
available both in symbolic computer algebra packages and as 
high-performance computing code, this opens up the possibility 
to calculate the luminosity density of the Sersic models to ar- 
bitrary precision. The wide range of analytical properties of the 
Meijer G function also allows to easily study the asymptotic be- 
haviour at small and large radii and to compute derivative quan- 
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tities such as the cumulative luminosity and the gravitational po- 
tential analytically. 

This paper is organized as follows. In Section 2 we discuss 
the general deprojection of the Sersic surface brightness profile 
using the Mellin transform method and we present a general ex- 
pression in terms of the Fox H function. In Section 3 we present 
more simple expressions for integer, half-integer and rational 
values of m in terms of the Meijer G function and discuss two 
special, interesting cases. In Section 4 we use the expressions 
for the luminosity density to calculate the total luminosity of the 
Sersic models, which serves as a sanity check on the derived 
formulae. In Section 5 we investigate the asymptotic behaviour 
of the luminosity density and in Section 6 we derive analytical 
expressions for the cumulative luminosity and the gravitational 
potential. Finally, in Appendix A we introduce the Meijer G and 
Fox H functions and discuss some of their most useful proper- 
ties. 



2. The luminosity density of the Sersic model as a 
Fox H function 

In spherical symmetry, the deprojected luminosity density v(r) 
at the spatial radius r, can be found from the surface brightness 
profile I(R) through the standard deprojection formula, 



v(r) 



i r 

n Jr 



dl 
d~R 



(R) 



dR 



If we introduce the reduced radial coordinate 

r 

S= K 



we find 



v(r) 



Rf, m 



X Js 



dt 



(2) 



(3) 



(4) 



The integral (4) cannot be evaluated in terms of elementary func- 
tions or even the standard special functions for general values of 
m. In order to evaluate it, we use a general method that builds on 
Mellin integral transforms and has become known as the Mellin 
transform method (Marichev 1983; Adamchik 1996; Fikioris 
2007). The Mellin transform 9JI/(m) of a function f(z) is defined 



as 



f(z)z'- 1 dz (5) 





The inverse Mellin transform is found as 



a«: 1 fe) = /(z) 



)z~" du 



(6) 



where the £, is a line integral over a vertical line in the complex 
plane. The driving force behind the Mellin transform method is 
the Mellin convolution theorem. The Mellin convolution of two 
functions /i(z) and /i(z) is defined as 

X°° It\ dt 

mh\j)- (7) 

Similar as for the well-known Fourier transform analogue, the 
Mellin convolution theorem states that the Mellin transform of a 
Mellin convolution is equal to the products of the Mellin trans- 
forms of the original functions, 



(8) 



Now it can be shown that any definite integral 

f(z)= g(t,z)dt 
Jo 



(9) 



can be written as the Mellin convolution of two functions f\ and 
fa. As a result, the definite integral (9) can be transformed to an 
inverse Mellin transform, 

f(z) = ^ fmfi(u)m fi _(u)z- u du (io) 

The power of this approach is that, if the functions f\ and fj 
are of the hypergeometric type, which is true for many elemen- 
tary functions and the majority of special functions, the inte- 
gral (10) turns out to be a Mellin-Barnes integral. Depending 
on the involved coefficients, this integral can be evaluated as a 
Fox H function, or in simpler cases, a Meijer G function (see 
Appendix A). 

The form of the equation (4) allows to apply the Mellin tran- 
form method, with z = 1 and 



R e mn 



and 



fi(t) = \ Vl - s 2 t 2 




if < t < s- 1 
else 



(11) 



(12) 



The Mellin transforms of these functions are readily calculated 

Iq 1 T(l +mu) 



R e 7T b" m 



l 



r(§) 



(13) 
(14) 



Substituting these values in the integral (10) and setting u = 2x, 
we obtain 



/o 



v(r) = - T s 



1 f i r T(l +2mx)T(\ + x) 



2ni 



X 



r(i +x) 



(b im s i y x dx 

(15) 

If we compare this expression with the definition (A. 11) of the 
Fox H function, we see that we can write the luminosity density 
of the Sersic models in compact form as 



f \ Io 1 -1 1/2,0 

V(r) = R e ^ S H ^ 



(1,1) 
(1,2m), (i 1) 



b 2m s z 



(16) 



3. Integer and rational values of the Sersic index 

Expression (16) represents a closed, analytical expression for the 
luminosity density of the Sersic model in terms of the Fox H 
function. While this function is receiving gradually more atten- 
tion both in mathematics and applied sciences, ranging from as- 
trophysics and earth sciences to statistics, its practical usefulness 
is still limited. In particular, no general numerical implementa- 
tions of the Fox H function are available (to our knowledge). 
Fortunately, the Fox H function can be reduced to a Meijer G 
function in many cases. In this Section we will derive an ana- 
lytical expression for the luminosity density of the Sersic model 
in terms of the Meijer G function for all rational values of m. 
The Meijer G function is much better documented and can be 
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an extremely useful tool for analytical work. The Meijer G func- 
tion has many general properties which allow to manipulate ex- 
pressions to equivalent forms, reduce the order for certain val- 
ues of the parameters, etc. Some of the most useful properties 
are listed in Appendix A, but there are many more. A particu- 
larly rich online source of information is the Wolfram Functions 
Site. 1 Moreover, commercial computer algebra systems such as 
Maple® and Mathematica® contain an implementation of the 
Meijer G function. It is also implemented in the open-source 
computer algebra package Sage and a freely available Python 
implementation of the Meijer G function to arbitrary precision 
is available from the Mpmath library. 2 

3. 1 . Integer and half-integer values of m 

If m is an integer or half-integer value, we can simplify expres- 
sion (16) by using the property 

N-i 

f]r(^) = ^-(2,)¥r ( z) (17) 

7=0 

Applying this recipe with N = 2m and z = 2m x gives 



r(l +2mx) = (2m)i +2mx (2ji) l 2-' n 



2m- 1 



Inserting this in expression (15), we find 



v(r) 



2Jo 

Re 



yfm _j 1 
(2ji) m S 2m 



i , \ 2m- 1 



2m I 



2m 



dx (19) 



Comparing this expression with the definition (A.l) of the 
Meijer G function, we obtain the following compact expression 
for the luminosity density of the Sersic model 



v(r) 



2Jo 

Re 



s G n n 

{2n) m °' 2m 



2m j 



2m 



with b a vector with 2m elements given by 

1 2 2m - 1 1 

'2 



b = 



2m 2m 



2m 



(20a) 



(20b) 



This expression is equivalent to (and actually even slightly sim- 
pler than) the expression obtained by Mazure & Capelato (2002) 
using the computer algebra package Mathematica®. Notice that 
Mazure & Capelato (2002) obtained their formula for integer 
values of m only, whereas our analysis shows that exactly the 
same expression also holds for half-integer values of m. 



1 The Wolfram Functions Site (http://functions.wolfram.com/) is a 
comprehensive online compendium, providing a huge collection of for- 
mulas and graphics about mathematical functions. It is created with 
Mathematica® and is developed and maintained by Wolfram Research 
with partial support from the National Science Foundation. A com- 
pendium of formulae on the Meijer G function can be found at 
http://functions.wolfram.com/PDF/MeijerG.pdf. 

2 Mpmath (http://code.google.eom/p/mpmath/) is a free pure-Python 
library for multiprecision floating-point arithmetic. It provides an exten- 
sive set of transcendental functions, unlimited exponent sizes, complex 
numbers, interval arithmetic, numerical integration and differentiation, 
root-finding, linear algebra, and much more. 



3.2. Rational values ofm 

Interestingly, these results can be generalized for all rational val- 
ues of m. Setting m = p/q with p and q integer numbers, we can 
write expression (15) as 



v(r) 



h 



1 

- f 

2m J £ 



qn\+2py)T{\+qy) 



r(l+qy) 



(b 2p s lq Y y dy (21) 



By multiple application of the identity (17) we can rewrite this 
expression in a format that leads to a Meijer G function. The 
result is 



y(r) = 



2/ ^JPI _j 2p+q-\,0 



Re 



q-l,2p+q-l 



2/' 



^1 



(2n)P 

with a a vector of dimension q - 1 with elements 

12 q-V 

q q'"" q 

and b a vector with 2p + q — 1 elements given by 
12 2p-l 1 3 2^-1 



b = 



2p ' 2p ' 



2p ' 2q' 2q 



2q 



(22a) 



(22b) 



(22c) 



It is straightforward to check that the general expression (22) 
reduces to the expression (20) for integer values of m. Using 
the order reduction formulae (A. 5) and (A. 6) of the Meijer G 
function, one can also demonstrate the expression (22) reduces 
to (20) for half-integer values of m. 



3.3. Special cases: m = 1 andm 



Among the family of Sersic models, there are two well-known 
specific cases for which the luminosity density can be calcu- 
lated explicitly in terms of elementary or special functions. The 
first of these two models is the exponential model, correspond- 
ing to m = 1. Exponential models are often used to describe 
the surface brightness profiles of dwarf elliptical galaxies (e.g. 
Faber & Lin 1983; Binggeli et al. 1984). If we introduce the no- 
tation h — R e /b, we can write the surface brightness profile as 



7(R)=/ e 



-R/h 



(23) 



If we deproject this surface brightness profile using the depro- 
jection formula (2), we recover the well-known result that the 
luminosity density can be written in terms of a modified Bessel 
function of the second kind, 



v(r) 



If we set m — 1 in the expression (20), we obtain 



v(r) 



J -l „2,0 

— s G, 
nR e 



0,2 





&V1 


1 1 

.2' 2 


4 



(24) 



(25) 



Using formula (A. 3), this expression reduces to the expres- 
sion (24). 

Another interesting special case is m = \, which cor- 
responds to a gaussian surface brightness profile. Such pro- 
files do not correspond to the observed surface brightness pro- 
files of galaxies, but they are very useful as components in 
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multi-gaussian expansions: even with a relatively modest set of 
gaussian components, realistic geometries can accurately be re- 
produced (e.g. Emsellem et al. 1994a,b; De Bruyne et al. 2001; 
Cappellari 2002). If we introduce <x = R e / y/2b and we use the 
total luminosity instead of the effective intensity as a parameter, 
we can write the surface brightness profile as 



I(R) 



(26) 



One of the key advantages of a multi-gaussian expansion of an 
observed surface brightness profile is that the corresponding lu- 
minosity density has a simple analytical form. Indeed, substitut- 
ing (26) into the deprojection formula (2), one can easily check 
that the deprojection of a gaussian distribution on the sky is a 
gaussian distribution as well, 



Kr)= (2^A eXp -i) 



(27) 



This result can also be found by setting m — | in equation (20), 



v(r) = s G 



0.1 



lbs 2 



(28) 



If we use equation (A. 2), we easily recover expression (27). 



4. The total luminosity 

The total luminosity of the Sersic model can be calculated by 
integrating the intensity on the plane of the sky, 



L = 2n 



r 

Jo 



I(R)RdR 



Inserting equation (1) one readily finds 



1 

b 2 " 



F(2m + 1) 



(29) 



(30) 



As a sanity check on the formula (22) and as a illustration of 
the power of the Meijer G function as an analytical tool, we can 
also calculate the luminosity by integrating the luminosity den- 
sity v(r) over the entire space, 



L-Aji\ v(r) r 2 dr 
Jo 

Inserting equation (22) we obtain 
, 2/ /?2 



(2n)P- 1 




-,2p+q-l,0 
T q-l,2p+q-l 



2 P 



2p 



(31) 



fi' 1 At (32) 



If we use the general property (A. 8) of the Meijer G function, 
we can evaluate this integral as 



2/0*2 [pi2 P ^n 2 j :r^{l^ J ) 



The product in the denominator can be simplified to 



(33) 



q-l 
7=1 



+ a 



iti) " r (i) 



(34) 



where the last transition follows from the identity (17). Similarly, 
the product in the numerator of equation (33) can be simplified 
to 



2p+q-l 



7=1 



^(2nr^-\2 P y- 2 i r(|) 



2q 



(35) 



If we substitute (34) and (35) into the expression (33), and we 
use m — p/q, we recover the expression (30) for the total lumi- 
nosity of the Sersic model, as required. 



5. Asymptotic behaviour 

One of the most useful properties of the general expression (22) 
is that it allows to elegantly determine the asymptotic behaviour 
of the luminosity density of the Sersic model at small and large 
radii. It is well-known (e.g. Ciotti 1991) that the Sersic models 
have a cusp for m > 1 and a finite luminosity density core at 
m < 1 ; this can be seen immediately by evaluating the integral 
(4) for r — (which converges only for m < 1), 



Vo 



= i±±- r 



e- b,i ""t^- 2 dt=— b'" T(l - m) (36) 

7lR e 



For a more detailed discussion on the behaviour of the lumi- 
nosity density at small radii, we can use the asymptotic expan- 
sion (A. 9) of the Meijer G function. In particular, this equation 
shows that the lowest order term in the expansion is determined 
by the smallest components bu in the vector b. This depends on 
the value of the Sersic index m. 



2q 



For m < l, the smallest component of the vector b is bi p 
and the second-smallest is b 



2 P +\ = Y q - After some algebra, 
which involves similar techniques as applied in Section 4, we 
find the asymptotic expansion 



v(r) 



hb_ 
ttR, 



r(i -m) i r(i 



-,l-3m 



3m) r. 
s 



(37a) 



For m — 1 the smallest component is still bi p = b% = \, but 
now the two components b\ and b^ p+ \ = b^ are both equal to 
j. In this case we cannot use the expansion formula (A. 9), since 
this formula is only valid if all components of the vector b are 
different. For m — | we find the expansion 



v(r) 

nR e 



2y + 21n/7-61n2) ^ 



(37b) 



where y w 0.57721566 is the Euler-Mascheroni constant. If \ < 
m < 1, bi p = y remains the smallest component of the vector 
b, but the second-smallest is now b\-j^. One finds 



y(r) 



hb_ 

7TR e 



F(l - m) y9r ^(2 2m) i_i 
b 1 -" 1 2m r(l 



j_ 

2ll! 



) 



(37c) 
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For m = 1, we again have two equal components in the vector 
b and we cannot readily apply formula (A. 9). The asymptotic 
expansion for small r now reads 



v(r) 



lob 



ln|-| - (y + \nb - In 2) 



(37d) 



Finally, if m > 1, the smallest component is Z7 L = which leads 
to 



v(r) 



* e 2»vsr(i-^) 



The five different asymptotic expansions (37) demonstrate the 
different behaviour of the luminosity density at small radii, de- 
pending on the value of the Sersic index m. For m < 1 the Sersic 
model has a finite density core with the central luminosity den- 
sity given by equation (36). At m — 1 the model has a logarith- 
mic cusp, at m > 1 we have a power-law cusp with logarithmic 
slope 1 — — . In particular, the de Vaucouleurs model has a lumi- 
nosity density profile that behaves as s~ 3/4 as small radii (Young 
1976; Mellier & Mathez 1987). Surprisingly, the Sersic models 
with m < 5 do not have a monotonically decreasing luminosity 
density profile with increasing radius. In the expansions (37a) 
and (37b), the first non-constant term has a positive coefficient 
and hence the luminosity density increases with increasing ra- 
dius in the nuclear region. The same accounts | < m < A, since 
the coefficient of the second term in the expansion (37c) is posi- 
tive for | < m < j and negative for A < m < 1 . 

At large radii, a single formula for the asymptotic expansion 
holds for all Sersic indices, 



, , lo b _ bs n„ I 1 

v(r) a / - — e 

R e V 2nm \ s 



in agreement with the result derived by Ciotti (1991). 



(38) 



6. Some other properties of the Sersic model 

The analytical expression (22) for the luminosity density of the 
Sersic models allows to express other properties of this family 
analytically in terms of the Meijer G function. The most im- 
portant ones are the cumulative luminosity and the gravitational 
potential. 

For a spherically symmetric system, they cumulative lumi- 
nosity L(r) can be calculated as 



L(r) = 4n 



I p(r') r' 2 dr' 
Jo 



(39) 



After substitution of expression (22) in equation (A. 7) we find 



L(r) = 



2hRl 



(2tt)p- 1 
which reduces to 



jP 2 Q2p+q-\,\ 

q q-Zp+q 



2p 



2p 



2I R 2 e 2 2 , 
L(r) = — — , — s G, 



(2n) 



m-1 



o 

b,-l 



\2m) ' 



(40) 



(41) 



for integer or half-integer values of the Sersic index m. 
Again, this expression is equivalent to the expression found by 
Mazure & Capelato (2002). The asymptotic expansion of the cu- 
mulative luminosity for small r can be found in the same way as 



we did for the luminosity density in Section 5. One finds after 
some calculation for m < 1 



L(r)~ -I Rib m r(l-m)s> 



For m — 1 we obtain 

AI R\b 
Uf) 



+ - y-lnfo + ln2 



(37e) and for m > 1 



Uf) 



2yfchR\b rQ-g;) 



.2+- 



2m + 1 



(42a) 



(42b) 



(42c) 



These asymptotic expressions can also be found by directly in- 
serting the equations (37) in formula (39). 

If we assume that mass follows light (or in case the Sersic 
model is used to describe the mass density), we can also calcu- 
late the (positive) gravitational potential ^(r). The most conve- 
nient way in the present case is to use the formula 



¥(r) = G Y 



f 



L(r')dr' 



(43) 



where the T is the mass-to-light ratio. The result reads 



{n ~ (2tt)p- 1 q^ S 



xG 



2[>+q,l 
q+\,2p+q+\ 



2q' q 



2p) 



(44) 



or if we introduce the total mass M — T L using expression (30) 
GM 



T(r) 



2p 

b 'i 



Re (2n)P^fpqY{ 2 f) 

1-i 1 



xG 



2p+q,\ 
q+\,2p+q+\ 



2q' q 



b 1 



(45) 



For integer and half-integer values of the Sersic index m, this 
expression simplifies to 



l2/» 



GM 



R e (2n) m aM r(2m) 



xG 



2m+l,l 
2,2m+2 



b \ 
2m ) 



(46) 



cient I appears in both the a and b coefficient vectors. Applying 



This expression can be reduced slightly further since the coeffi 
l 

2 

equation (A. 5), the final result reads 
GM b 2m 



sG 



2m, 1 
l,2m+l 



R e (2n) m yfm F(2m) 
with V a vector with 2m + 1 elements given by 
1 2 2m - 1 1 



2m 



2m 



b' 



2m ' 2m 



2m 



(47a) 



(47b) 



This expression is somewhat simpler than, but equivalent to ex- 
pression (28) in Mazure & Capelato (2002). 
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Since the luminosity density of the Sersic models never falls 
steeper than r _] at small radii, it is no surprise that all Sersic 
models have a finite potential well for all values of to. Ciotti 
(1991) already derived an expression for the depth of the poten- 
tial well using the general expression 



¥ = -AG T 



r°° ai 



dR 



(48) 



Applied to the Sersic model surface brightness profile, the result 
reads [his equation (12)] 



V = AGri u R t 



F(to + 1) 



Taking the limit r 



for the expression (45), we find 
GM 2b m T{m) 



^0 = 



R e 7T T(2m) 



(49) 



(50) 



equivalent to (49). Using the expansion formulae for the Meijer 
G function, we can calculate the asymptotic expansion for the 
potential at small radii. Not surprisingly, we again obtain differ- 
ent expansions for m smaller than, equal to, and larger than 1. 
After a lengthly calculation, one finds for to < 1 a quadratically 
decreasing potential, 



<F(r) 



GM 



2b m F(to) 



'3m 



r(l-m) 2 
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7i T(2m) 3mn Y(2m) 
For the exponential model m = 1, one obtains 



(51) 
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(52) 



where i[/ z is the digamma function. For to > 1 the potential de- 
creases more softly than quadratically, 



Y(r) 



GM 

~R~7 



2b m Y(m) 



n T(2m) 
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J+ 1 



(53) 



Finally, at large radii, the potential of all Sersic models falls off 

as 

<P(r) ~ — (54) 
r 

as required for a system with a finite mass. 



7. Conclusions 

We have used the Mellin transform technique to derive a closed, 
analytical expression for the spatial luminosity density v(r) of 
the Sersic model. For general values of the Sersic parame- 
ter to, this expression is a Fox H function. We derive simpli- 
fied expressions for v(r) in terms of the Meijer G function for 
all rational values of to; for integer values of to our results 
are equivalent to the expressions found by Mazure & Capelato 
(2002). Our analytical calculations complement other theoreti- 
cal studies of the Sersic model (Ciotti 1991; Ciotti & Lanzoni 
1997; Ciotti & Bertin 1999; Trujillo et al. 2001; Cardone 2004; 
Graham & Driver 2005; Elfasdottir & Moller 2007) and, given 



the extended literature on the analytical properties of the Meijer 
G function, can be used to further examine the properties of 
this model analytically. We have investigated the asymptotic be- 
haviour of the luminosity density at small and large radii, and 
find a rich variety in behviour depending on the value of to. We 
have also derived analytical expression for derived quantities, 
in particular the cumulative luminosity and the gravitational po- 
tential. Our results can also be used in practical calculations: as 
implementations of the Meijer G function are nowadays avail- 
able both in symbolic computer algebra packages and as high- 
performance computing code, our results open up the possibility 
to calculate the luminosity density of the Sersic models to arbi- 
trary precision. 
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Appendix A: The Meijer G and Fox H functions 

The Meijer G function is a universal, analytical function that has been introduced as a generalization of the hypergeometric series 
by Meijer (1936). Nowadays it is more commonly defined as an inverse Mellin transform, i.e. a path integral in the complex plane, 



/~*m,n 
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b it . 



1 r in: 

~2mJ £ U% m+l 
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z~ s ds 



(A.l) 



with _£ a path in the complex plane and a and b are two vectors of dimension p and q respectively. 

The general Meijer G function can reproduce many commonly used special functions, including Bessel functions, elliptic inte- 
grals and hypergeometric functions. Some of the special cases are 
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(A.2) 
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where K v {z) is the modified Bessel function of the second kind. 

The Meijer G function has numerous useful properties that allow to transform expressions to equivalent expressions. For exam- 
ple, one can shift all parameters by a given number, 
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An important property is that one can reduce the order of the Meijer G function if one parameter appears in both the upper and lower 
parameter vectors (depending on the position). For example, if one of the a k with n < k < p equals one of the bj with 1 < j < m, 
then 
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Another powerful property that allows to reduce the order of the Meijer G function in certain cases is the following 
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where k is a positive integer number, 6 and v are defined as 

p + q 
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(A.6b) 



and the vectors a' and b' are defined as 
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One of the most powerful properties of the Meijer G function as an analytical tool is that several integrals involving Meijer functions 
can be evaluated in terms of higher-order Meijer functions. For example, 
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The corresponding definite integral can be evaluated as follows 
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Another useful property is the asymptotic expansion of the Meijer function for small z, in the case of p < q and simple poles, 



°P,q 



a 




b 


z 



wu^nbj-h) n ^r(i-fl J+ ^) h 

naj-h) U% m+ J(l-b j + b k ) Z 



n 



k=\ 1 i 7=n+l 



I _|_ f_y\-m-n+p 



U P j=i a-aj + b k ) 
U% l (l-b J + b k ) 



z + 0(z 2 ) 



(A.9) 



For the asymptotic expansion at large z, one can use the identity 
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The Fox H function (Fox 1961) is a generalization of Meijer G function and is also defined as an inverse Mellin transform, 
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Not surprisingly, the Fox H function shares many of the properties of the Meijer G functions and complete volumes have been written 
about its identities, asymptotic properties, expansion formule and integral transforms (e.g. Mathai & Saxena 1978; Kilbas & Saigo 
2004; Mathai et al. 2009). 



